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axisymmetric systems 
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ABSTRACT 

In a previous paper the complex-shift method has been apphed to self-gravitating 
spherical systems, producing new analytical axisymmetric density-potential pairs. We 
now extend the treatement to the Miyamoto-Nagai disc and to the Binney logarithmic 
halo, and we study the resulting axisymmetric and triaxial analytical density-potential 
pairs; we also show how to obtain the surface density of shifted systems from the 
complex-shift of the surface density of the parent model. In particular, the systems 
obtained from Miyamoto-Nagai discs can be used to describe disc galaxies with a 
peanut-shaped bulge or with a central triaxial bar, depending on the direction of the 
shift vector. By using a constructive method that can be applied to generic axisym- 
metric systems, we finally show that the Miyamoto-Nagai and the Satoh discs, and the 
Binney logarithmic halo, cannot be obtained from the complex-shift of any spherical 
parent distribution. As a by-product of this study we also found two new generating 
functions in closed form for even and odd Legendre polynomials, respectively. 
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1 INTRODUCTION 

For the discussion of many astrophysical problems where 
gravity is important, a major difficulty is set by the poten- 
tial theory, in particular in the applications where the avail- 
ability of analytical density-potential pairs deviating from 
spherical symmetry is important. Due both to its relevance 
and its relation with electrostatic and magnetostatic (see 
Jackson 1999), it is natural that different techniques have 
been developed to study this problem (e.g., see Binney & 
Tremaine 2008; see also Ciotti & Berlin 2005; Ciotti & Gi- 
ampieri 2007, hereafter CG07, and references therein). 

In this paper we focus on a particular but very elegant 
technique. In fact, in the context of classical electrodynamics 
Newman (1973; see also Newman & Janis 1965, Newman et 
al. 1965) considered the case of the electromagnetic field of 
a point charge displaced on the imaginary axis; the remark- 
able properties of the resulting field have been studied by 
Carter (1968), Lynden-Bell (2000, 2002, 2004a,b, and refer- 
ences therein; see also Teukolsky 1973; Chandrasekhar 1976; 
Page 1976), and Kaiser (2004, and references therein). This 
procedure, also known as the complex-shift method, was first 
introduced by Appell (1887; see also Whittaker & Watson 
1950) to the case of the gravitational field of a point mass, 
and successively used in general relativity (e.g., see Gleiser 
& PuUin 1989; Letelier & Oliveira 1987, 1998; DAfonseca, 
Letelier & Oliveira 2005 and references therein). In the case 



of classical gravitation, this method has been extended to 
continuous mass density distributions by CG07, who consid- 
ered spherical parent systems, and obtained new analytical 
axisymmetric density-potential pairs. 



This paper is the natural follow-up of the study of CG07 
to axisym-metric parent systems, not only to further extend 
the class of new axisymmetric systems, but also to obtain 
new triaxial density-potential pairs. In Section 2 we briefly 
recall the framework developed in CG07, and then we dis- 
cuss the application of the method to axisymmetric parent 
models. In Section 3 we describe the new axisymmetric and 
triaxial systems obtained from the axial and the equatorial 
shift of the Miyamoto-Nagai disc and the Binney logarith- 
mic halo, while in Section 4 we discuss the projected surface 
density of the new models. Section 5 summarizes the main 
results obtained. In Appendix A we discuss the expansion 
of the complexified potential of a generic assigned density 
distribution in terms of the shift amplitude and, as a by- 
product, we found two new generating functions in closed 
form for even and odd Legendre polynomials, respectively. 
In Appendix B we finally present the general solution of an 
inverse problem associated with the complex-shift method. 
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2 GENERAL CONSIDERATIONS 

Here we recall the basic idea of CG07, who extended the 
complexification of the Coulomb field of a point charge (e.g., 
see Lynden-Bell 2004b), to the gravitational potential $(a;) 
of a (nowhere negative) density distribution p{x). From 
now on a; = (a;, y, z) indicates the position vector, while 
< x,y >= XiXji is the standard inner product over the reals 
(repeated index summation convention implied) . Let assume 
that p{x) satisfies the Poisson equation 

V^-I- = 4-!vGp, 

and that the associated complexified potential f&c with shift 
la is defined as 

^c{x) = ^{x — m), 

where i^ = — 1 is the imaginary unit and a is a real vector. 
From the linearity of the coordinate transformation and of 
the Laplace operator it follows that 



V^'^c 



where 



iTvGpc 



pc{x) = p{x ~ la) 



is the shifted density 

|2 ll_l|2 



we recall that in this technique 
x — ia\\~ = ||a^ir — ||a|P — 2i < a,x > (e.g., see Lynden- 
Bell 2004b). By separating the real and imaginary parts of 
$c and pc obtained from the shift of a known real density- 
potential pair (or, if simpler, taking the Laplacian of $c), one 
obtains two new density-potential pairs. Of course, positivity 
of the new densities is not assured. In fact CG07 proved that 
for a generic parent density p, the imaginary part Q'(pc) nec- 
essarily changes sign because the total mass (and the total 
gravitational energy) of the complexified density-potential 
pair coincides with the corresponding quantities of the par- 
ent (real) density-potential pair, so that (from the mass in- 
variance) it follows that J 9(pc) <jPx — 0. This implies that 
Sy(pc) cannot be used to describe a gravitating system. On 
the other hand, either cases can happen for the real part 
5R(pc), depending on the specific seed density and shift vec- 
tor adopted. For example, while iR(pc) of the shifted Plum- 
mer (1911) and Isochrone (Henon 1959) spheres is positive 
over the whole space when the module of the shift vector is 
less than a critical value Umax, 5R(pc) of the shifted singular 
isothermal sphere has negative regions independently of the 
adopted a 7^ (see CG07). 



2.1 Axisymmetric parent systems 

At variance of the spherical case, the result of the com- 
plex shift of an axially symmetric potential ^{R,z), where 
R = \/x^ + y^ is the cilindrical radius, depends on the di- 
rection of the shift vector a. For simplicity in the follow- 
ing Sections we consider two special shift directions, namely 
the axial shift a = (0,0, a), and the equatorial shift, where 
without loss of generality we assume a = (aa;,0, 0); how- 
ever, some general results for arbitrary shift direction and 
parent density distributions are postponed to Section 4 and 
Appendix A. The axially shifted potential and density are 



(2) 
(3) 
(4) 



where 

2 ; 

2(. = 2 

while in the case of equatorial shift 

$c = ^(-RcZ), Pc=p(-Rc,2), 

where 

R^ — R — a — 2iax 

is the shifted cylindrical radius. 

3 RESULTS 

3.1 The shifted Miyamoto-Nagai disc 

We start by considering the family of models originated by 
the complex shift of the Miyamoto-Nagai disc (1975, here- 
after MN; see also Binney & Tremanine 2008). For conve- 
nience we use the relative potential 'I' = — "I> so that the MN 
density-potential pair can be written as 

(5) 



^{R,z)-- 


VR^ + {s + 0^ C ^ 


p{R,z) = 


/ 7s 5s2 si?2 3 '^ 



(6) 



In the expressions above ^ = y'T+z^, all the coordinates 
are in units of the scale-length b, while the relative poten- 
tial and the density are normalized to GM/b and to M/b'^, 
respectively; M is the total mass of the system. For increas- 
ing s the density distribution becomes more and more flat 
and concentrated, while for s = it reduces to the Plummer 
sphere, for which CG07 has already discussed the complex- 
iflcation and showed that for a < 0.588 the real part of the 
shifted density is nowhere negative. We then expect that 
the method applied to MN discs will again produce nowhere 
negative real densities, at least for discs with low flattening. 

3.1.1 The axial shift 

The axially shifted MN density can be obtained by the direct 
complexification of eq. ([6]) as 



Pc 



o 7s 5s2 



sR^ + s' 



47r ' 



(7) 



in the following the shift modulus a is intended expressed 
in units of b. We start from the evaluation of ^c. From ^^ = 



V(l-a2-h22)2_^4„222 



sin^ : 



2a2 
u 



(8) 



(9) 



where in eq. ^ (and in all analogous cases) the square root 
is the positive arithmetic operator; thus 



Cc 



9/2+ikT: 



(fc = 0,l). 



(10) 



$c=$(i?,2c) 



Pc = p{R,Zc), 



(1) 



Note that cos 9 > for a < 1: for simplicity we restrict 
to this case. In particular, we cut the complex (u, 9) plane 
along the negative real axis so that —n < 8 < tt.Iu addition. 
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Figure 1. Top panels: 5ft(pc) of the axially shifted MN disc with s = 3/2. Left and central panels; logarithmically spaced isodensity 
contours on a meridional plane of the parent (a = 0) and of the shifted model with a = 3/4, respectively. Top right panel: edge-on surface 
density for the a = 3/4 model. Bottom panels: 3f?{pc) of the axially shifted Binney logarithmic halo with q = 0.71. Left and central 
panels: logarithmically spaced isodensity contour on a meridional plane of the parent (a = 0) and of the shifted model with a = 1/2, 
respectively. Bottom right panel: edge-on surface density for the a = 1/2 model. 



the natural request that ^c > for z = leads to choose 
fc = 0. The bisection formulae for 6 are 



cos - — 
2 



1 + cos 61 . 9 _ sin 6* 

2 ' '*''' 2 " 2cos(6l/2) 



(11) 



where the sign determination is dictated by the interval of 
variation of 6, and the real and imaginary parts of ^c are 
obtained from eqs. (|8))- (|1H) as 



mc) = 



1 — a^ + z'^ + u 



3(?c) = - 



mc 



(12) 



We are now in position to evaluate $c = 1/Cc from eq. 
((5l). By writing t^^ = R^ + {s + Cc)^ = ve"^, we obtain after 
some algebra 



R-" 



ucos6 + 2s\/ucos ■ 



sintp = 



usin0-|-2s^sin(6i/2) 



(15) 



From the positivity of cos 6 it follows that also costp > 0, 
while sin;^ changes sign as sin6': the formulae for cos{(p/2) 
and sm{(fi/2) are analogous to those in eq. (Ill|l . and the 
same considerations followed to determine ^c apply to (c, so 
that (c = %/u6"^ and 



*c 



"fi/2 



Co 



(16) 



The real and imaginary part of the shifted potential are 
given by 



5R(*c) = 



^l - a'^ + z^ + V + K-^ + s'-^ + 2s3ff(Cc) 



(17) 



-|- [usind + 2s\/usin- 



1/2 



COSip — 



R^ + s'^ +UCOS9 + 2a^cos{e/2) 



(13) 



(14) 



9(*c) = 



vm{-^^ 



m^ 



(18) 



The complexified density can now be easily determined: the 
simplest way is to substitute ^c (eq. [10] with k = Q) and ^c 
(eq. [16]) in eq. ([7j), and to use eqs. ([T3j)-([l5l) and lfr7|l - (fT8)) 
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Figure 2. Top panels: logarithmically spaced isodensity contours of 5R(pc) of tlie equatorially sliifted MN disc with s = 3/2 and 
a = 67/40, in the three orthogonal coordinate planes. The shift is along the x-axis and all the lengths are in units of b. Central panels: 
logarithmically spaced isodensity contours of SR(pc) of the equatorially shifted Binney logarithmic halo with q = 0.71 and a = 0.85, in 
the three orthogonal coordinate planes. The shift is along the x-axis. Bottom panels: the corresponding projected surface densities of the 
Binney logarithmic halo. 



after separation. The explicit formulae are however quite 
cumbersome and are not reported here. In any case, a first 
check was obtained by comparing term by term the limiting 
case s = with the formulae derived in CG07 for the shifted 
Plummer sphere. As an additional check, we also compared 
the Laplacian of 3fJ(^c) and Q'(*I'c) with the corresponding 
components of the shifted density, for generic s. 

As expected, Q'(pc) changes sign crossing the equatorial 



plane of the system. For what concerns 3fJ(pc), a numerical 
exploration reveals that the density becomes negative on 
the symmetry axis at 2 ~ 0.811 for s — 1/2 and a > 0.672, 
while the same thing happens at z ~ 0.742 for s — 3/2 
and a > 0.764. Thus, Umax increases for increasing flatten- 
ing. In Fig. [T] (first and second panels of the top row) we 
show the isodensity contours of 3fJ(pc) in a meridional plane 
for the unshifted model and the model with a — 3/4, i.e. a 
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value of the shift parameter near to the critical value. It is 
apparent how at large distance from the centre the density 
distribution is almost unaffected by the shift: this is not a 
pecularity of the present model, but it is a general prop- 
erty of the complex shift, as already recognized by CG07 
in the spherical cases, and proved in general in Appendix 
A. The second feature is the toroidal shape of the density 
near the centre, with a characteristic depression along the 
shift axis: such structure is curiously similar to the toroidal 
shapes of the homeoidally expanded density-potential pairs 
of Ciotti & Bertin (2005) and to MOND density-potential 
pairs (Ciotti, Nipoti & Londrillo 2006), obtained through 
potential mapping of homeoidally expanded Newtonian dis- 
tribution. The resulting density distribution is reminescent 
of peanut-shaped bulges (e.g., Shaw, Wilkinson & Carter 
1993; Kuijken & Merryfield 1995; Emsellem & Arsenauh 
1997; Bureau & Freeman 1999; Liitticke, Dettmar & Pohlen 
2000), and the axially shifted MN disc could then be used 
as a compact analytical density-potential pair to investigate 
stellar orbits and gas dynamis in disc galaxies with peanut- 
shaped bulges (e.g., see Binney & Petrou 1985; Combes et 
al. 1990; Patsis, Skokos & Athanassoula 2002; Quillen 2002; 
Athanassoula 2005; De Battista et. al 2005). 

3.1.2 The equatorial shift 

We now consider the equatorial shift of the MN disc and 
without loss of generality we assume the shift direction along 
the X axis. In this case the variable ^ in eq. ([5]) is unaffected 
by the shift, while the variable R transforms according to 



°, where 



(19) 



(20) 



eq. ^ as Cc = R^ + is + (.f - a^ - 2\ax 
v^ x/[i?2 + {s + ^)2 _ a2]2 + 4a2a;2 , 



H^ + (s + n^ - a^ . 2ax 
cos (p = , sm ip — , 

V V 

and 

C, = V^e'^/'+"'" , (fc = 0, 1) . (21) 

For simplicity we restrict to the case a < 1 + s, so that 
cosip > everywhere. Again, we cut the complex {v,ip) 
plane along the negative real axis, and we fix fc = so that 
Cc > for z = 0. Accordingly, cos(iy9/2) and sin(93/2) are 
given by identities analogous to eq. (|11|) . and simple algebra 
shows that 

y^v + R2 + {s + 0^ - a2 



K(*c) 



9(*c) = 



vV2 



(22) 



(23) 



«2K(*c) 

As for the axial shift, also in the case of equatorial shift the 
the simplest way to obtain the shifted density is by direct 
complexification and expansion of eq. (|6]). In particular, the 
function ^ is unaffected by the shift, and after some work 
we found 



^{Pc) = 



o 7s 5s^ 

3+ J + ^ + 



s{R'' + s^ 



An 



where 



3fJ(*c 



and 



Sjf*n = 



3fJ(*c 



55R(*c 



lOa^a;^ 




Sa-'x* 1 


v^ 


+ 


U83fj(*c)^ 


lOa'^x 


2 


a^x* 



3ff(*c 



«8K(*c)* 



K*c). 



(25) 



(26) 



The obtained formulae have been verified by evaluating an- 
alytically the Laplacian of eq. (|22p and comparing the re- 
sulting formulae with eq. (|24p . In the special limit s = the 
new pair coincides, as expected, with the formulae of CC07 
for the shifted Plummer sphere. 

For what concerns 5R(pc), numerical explorations reveal 
that for s = 1/2 the density 5R(pc) becomes negative on the 
X axis at a; ~ 1.069 and a > 0.943, and at a; ~ 1.563 and 
a> 1.699 for s = 3/2. The values of amax and s are corre- 
lated as in the case of axial shift, in the sense that flatter 
models are able to support larger shifts. In the three top 
panels of Fig. [51 we present the isodensity contours of 5R(pc) 
for s — 3/2 and a near its critical value in the three or- 
thogonal coordinate planes. The effect of the shift along the 
X direction is apparent. In particular, note that isodense in 
the top left panel (equatorial plane section) would be per- 
fectly circular for the parent unshifted MN disc. The same 
striction effect is apparent in the top right panel (section on 
the meridional plane containing the shift vector). Finally, 
the top central panel is the meridional section in the plane 
orthogonal to the shift vector and the density reminds that 
of the parent MN disc. In fact, it is easy to prove that in 
this plane the isodense corresponds to that of a standard 
MN disc where y^ is replaced by y^ — a^. Thus, the resulting 
model is a new triaxial density-potential pair, with a triaxial 
central body smoothly joined to an axisymmetric disc in the 
external regions. Models of this family, due to their simple 
analytical nature, could be used to study how central bars 
influence the orbits of stars and gas in galaxy discs, and 
to provide the gravitational fleld needed in hydrodynami- 
cal simulations of black hole accretion at the centre of disc 
galaxies via bar fueling. 

3.2 The shifted Binney logarithmic halo 

We now focus on the Binney (1981) logarithmic halo (see 
also Binney & Tremaine 2008) 

Infl + m^) , , ^ 

y I ' -' " (27) 



$ = 



2 _ r>2 , ^ 

m = R ~\ — - 



2 ' g- 

2 [1 + R\l - q^)] e-"" + (2g^ - l)e-^^ 
47rg2 



2[q^ + z\l-q-')]e- 



+ e- 



47rg2 



(28) 



where all the coordinates are normalized to a scale-lenght 
b, the potential to vq (the asymptotic equatorial circular 
velocity), and the density to vl/Gh^ (Evans 19930- The 
parameter q > l/\/2 controls the model flattening, and the 



as2;9(\['c) 



(24) ^ Note that in oq. (23) of Evans (1993) the quantity R^ is missing 

in the last term and i?7 in the dimensional coefficient. 
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lower limit is required by the positivity of p. In the following, 
we will focus on the real part of the shifted pair, as we 
already know that the pair Q{pc)-Q{^c) cannot be used to 
describe a gravitating system. 

3.2.1 The axial shift 

From eq. ((2} it follows that the axially shifted potential is 



-l>c = iln(^l^ 


2 

a 


so that 




K($c) = i In 


[('- 



2ia2:\ 



>-4 



+ 



ia^z^ 



(29) 



(30) 



The real part of the complexified density can be obtained 
from the Laplacian of eq. (I30|l , however we use here the first 
identity in eq. (I28|) . In fact, as the variable R is unaffected by 
the axial shift, we just substitute eq. (|29p in the expression 
for p{R^ , ^) obtaining after some work 



^(Pc) = 



+ 



(2g2 



1] 



47rg2 



,4sre(Sc) 



40^22 



-4!R(Sc) 



27rg2 



,4»{4.c) 



-SSRC-Sc) 



(31) 



As in the other cases, we verified by Laplacian evaluation 
the correctness of eq. (|3Hl . Note that solving eq. (|30p for 
the variable 2^, eq. (I31|l can be recast in an algebraic form 
5R(pc) = F[R ,5R("I>c)], i.e. in a form suitable (in principle) 
for the recovering of the even part of the two- integrals phase- 
space distribution function (e.g., see Fricke 1953; Lynden- 
BeU 1962; Toomre 1982; Hunter & Quian 1993; Ciotti & 
Berlin 2005). Unfortunately, the presence of (algebraic) ir- 
rationalities seems to exclude the possibility of a simple an- 
alytical inversion. 

As for the axial shift of MN discs, the global positivity 
of 5R(pc) cannot be guaranteed for a generic value of a. In- 
deed, numerical explorations show that for g = 1 (the parent 
spherical case) the density becomes negative on the simme- 
try axis of the system at 2 ~ 0.464 for a > 0.808, and for 
g = 0.71 (the parent maximum fiattening) a negative region 
appears at 2 ~ 0.499 for a > 0.503. Note that at variance 
with the MN case, the critical a decreases for parent sys- 
tems increasingly deviating from spherical symmetry (i.e., 
with smaller g). This is not surprising: in fact, while for MN 
models the flatter models are more dense on the symmetry 
axis, the logarithmic halo develops a toroidal shape, which is 
less able to support the additional density decrease produced 
by the shift. In Fig. [T] (bottom panels) we show the isoden- 
sity contours in a meridional plane of K(pc) for the model 
with a shift parameter near the maximum: the appearence 
of the critical region on the symmetry axis is evident. 

3.2.2 The equatorial shift 

The equatorial shift of potential (|27|) is obtained from eq. 
0, leading to 



•l-c^ 



In (1 - a^ + m^ - 2iax) 



(32) 



and 



3fJ($c) 



[n[{l-a^ + Tn^y+4a^x^] 



(33) 



Again, the real part of the complexified density can be ob- 
tained by direct evaluation of the Laplacian of eq. (I33|l , but 
here we use the second identity in eq. (|28p , so that only the 
evaluation of p(z^, $c) is needed. The real part of complex- 
ified density is then given by 



5R(pc 



(l-a^+m^ 



47rg2 



+ 



[q' + ^'{l 



-2)][e4K(*.) 



27rg2 



(34) 



The positivity of 5ft(pc) cannot be guaranteed every- 
where for a generic value of a because a negative term is 
present in its expression. Numerical explorations show that 
this term dominates over the positive for g = 1 at 2; ~ 0.464 
when a > 0.808, being this case coincident with that of the 
axial shift of the spherical model, and for g — 0.71 at 
X ~ 0.372 when a > 0.853. Note that, in contrast with the 
axial shift case, Umax decreases for increasing q. This is ex- 
pected from the presence of density lobes on the equatorial 
plane of the parent distribution. In Fig. [2] (central panels) we 
present the isodensity contours of 5R(pc) in the case g = 0.71 
and a = 0.85, in the three orthogonal coordinate planes. The 
triaxial nature of the resulting structure is apparent, as the 
similarity of the density in the central panel with that of the 
parent system, already mentioned in the MN case. 



4 SURFACE DENSITIES 

So far, we have shown how the complex-shift method can be 
applied to the case of axisymmetric systems to obtain new 
analytical density-potential pairs. In this Section we extend 
the analysis to the projected surface density of the shifted 
models. In general, the surface density Ec of a complexified 
system projected along the direction of n (unit vector) is 
defined as 



Ec(a 



/+00 
-00 



Pc{xj_ +ln)dl , 



(35) 



where I —< x,n > is the line-of-sight coordinate, and 
x± = X — In is the position vector in the projection plane 
< x±,n >— 0. From the linearity of the projection inte- 
gral, it follows that the real and the imaginary parts of Ec 
are the projections of the real and imaginary parts of the 
shifted density. Remarkably, Ec can be obtained directly by 
shifting the surface density E of the parent system with 
shift a± — a— < a, n > n. In fact, pc ~ p{x — la) = 
p[x± — ia± + {l — i < a,n >)n], and the integration along I 
in eq. (|35|) is unaffected by the change of the origin, as the 
integration interval extends from —00 to +00. Thus, from 
eq. (|35|l it follows that 



Ec(a 



= E(xx — lax) 



(36) 



in particular, the surface density of a model projected along 
the shift direction {a± = 0) is real and coincides with the 



surface density of the parent density distribution. For sim- 
plicity, in this Section we will limit our attention to the face- 
on and the edge-on projection of the models constructed in 
Sect. 3. 



4.1 Edge-on projection 

The edge-on projection of an axisymmetric density p{R, z) 
can be written without loss of generality as 



E(a;, 



p{\/x^ + y^,z)dy. 



(37) 



In the present context we identify two special configurations, 
i.e. the edge-on projection of an axially shifted system, and 
of a system shifted equatorially in the direction perpendic- 
ular to the line-of-sight so that, according to eq. (|36p . 



Ec(a;, z) — T,(x, z — la) 



Ec(a;, z) = E(a; — la, z) 



(38) 



respectively. 

As a first example we consider the edge-on projection 
of the MN disc, for which the surface density (normalized to 
M/6^) can be written as 



T,(x,z) = 



where 



A* 



2 H 1 1 

^ e ?^ / 27r 



(39) 
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4.2 Face-on projection 

The face-on projected surface density of an axisymmetric 
density p{R,z) is given by 



(42) 



E(_R) = / p{R,z)dz, 

so that in the case of axial shift 

Ec(i?) = E(i?) , 

while in the equatorial shift a = a±, and 

Ec(i?) = E(i?c) . 



(43) 



(44) 



As an example we show the face-on surface density of 
the equatorially shifted Binney logarithmic halo (normalized 
to vf/Gb) 

g(2 + i?2) 



E{R) 



(45) 



4(l+it:2)3/2 ' 

where the radius is in units of b and the complexification is 
obtained from eq. Q. Isodensity contours are represented 
in Fig [5] (bottom left panel), for a system with a flattening 
parameter q — 0.71, and shift along the a;-axis of a = 0.85. 
It is apparent how the spatial deformation is mitigated by 
the projection. Unfortunately, the face-on projection of the 
MN disc cannot be expressed in closed form, so that we do 
not discuss this case. 



A = 



^x^ + {s + 0' 



(40) 



(Ciotti & Pellegrini 1996). From the general considerations 
above, it follows that eq. p9|l is also the edge-on surface 
density of a MN disc equatorially shifted along the line-of- 
sight direction. In the case of axial shift, the complexification 
is obtained by shifting ^ — > ^c and A — > Ac in eq. (I39p . 
according to the first identity in eq. p8p . Note that the 
function A is the potential evaluated on the projection plane 
(y = 0), so that one can use in eq. (|39p the quantities ^c and 
^c in their polar form (eqs. [10) and [TB], respectively). In 
the case of shift along x, the function ^ is unaffected, while 
Ac is obtained as the equatorially shifted potential (eq. [22] ) 
evaluated on the projection plane y = 0. 

For simplicity in Fig. [l] (top right panel) we show only 
the edge-on case of the axial shift of MN disc. The peanut 
shaped contours are clearly visible, even though less pro- 
nounced than in the spatial section (top central panel), as 
the natural consequence of projection. 

The edge-on projection of the Binney halo is elementary, 
and the surface density, normalized to v^/Gb, is 



T,{x,z) = 



-x^ + z^ + q'^ 



4g2 (1 + a;2 + z2/g2)3/2 



(41) 



also the complexification is trivial, and the resulting formu- 
lae are not reported here. The edge-on view of the axial 
shift is given in Fig. [l] (bottom right panel) , while in Fig. 
[5] (bottom central and right panels) we show the edge-on 
projection for a equatorially shifted model along the line-of- 
sight (central panel), and perpendicular to it (right panel). 
Again, the general considerations of Sect. 4 assure that the 
surface density in the central panel coincides with that of 
the parent unshifted halo. 



5 CONCLUSIONS 

In this paper we have shown that the complex-shift method 
is able to produce new and explicit density-potential pairs 
with finite deviation from spherical symmetry. In particular, 
the imaginary part of the complexified density corresponds 
to a system of null total mass, while the real component can 
be positive everywhere (depending on the original parent 
density distribution and the amount of the complex shift), 
and so used to describe astrophysical systems. 

In a natural follow-up of the preliminary investigation 
of CG07 (that was restricted to parent spherical models) 
here we considered the case of axial and equatorial shift 
of axisymmetric systems. As illustrative examples we con- 
sidered the cases of the Miyamoto-Nagai disc and Binney 
logarithmic halo, and we found that the shift leads to new 
and explicit density-potential pairs, with everywhere posi- 
tive density (for a shift vector modulus smaller than some 
limit value). From a qualitative point of view, the shifted 
densities are more and more similar to the parent distri- 
bution at large and larger distances from the centre, and 
this behaviuor is elucidated in Appendix A by means of the 
integral representation of the complexified potential. As a 
by-product of this approach we found two generating func- 
tions in closed form for even and odd Legendre polynomials, 
respectively, that at the best of our knowledge were not pre- 
viously known. 

We also showed how the projected surface density of a 
generic complexified density can be obtained by shifting the 
projected density of the parent distribution, thus avoiding 
the projection of the new density. This would be in general 
more difficult, due to loss of symmetry as a consequence of 
the shift action. As an application, we constructed the sur- 
face densities of the new models presented. In particular. 
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we found that the edge-on surface densities of the axially 
shifted Miyamoto-Nagai disc resemble the surface brightness 
of the so-called peanut-shaped bulges. Also relevant is the 
case of equatorial shift, leading to a family of simple analyt- 
ical density-potential pairs axisymmetric at large distances 
but triaxial in the central regions. Thus, the systems pro- 
duced by the complex shift of Miyamoto-Nagai discs, could 
be used as simple analytical models to study gas and orbital 
dynamics in axisymmetric peanut-shaped galaxies. 

In this paper, for simplicity, we restricted the anal- 
ysis to some simple axisymmetric density-potential pairs, 
but the class of the systems obtained via the complex-shift 
method could be extended by complexification of the gen- 
eralized Miyamoto-Nagai and Satoh density-potential pairs 
constructed by Vogt & Letelier (2005, 2007) or by a joint 
application of the complex-shift and Kelvin inversion to the 
flat rings model of Letelier (2007). 

We conclude by discussing an inverse problem related 
to the present investigation. In fact, as shown in CG07, the 
real part of a complex-shifted spherical system corresponds 
to a density-potential pair characterized by axial symmetry. 
In the course of an interesting discussion, Tim de Zeeuw 
suggested to investigate wheter the most common axisym- 
metric density-potential pairs in the literature could be the 
real part of some complexified spherical system. Quite sur- 
prisingly, it turns out that it is possible to formulate a sim- 
ple and constructive approach (presented in Appendix B) 
able to solve this question in general. In our cases, we were 
able to prove that the three investigated models, namely the 
Miyamoto-Nagai and Satoh discs, and the Binney logarith- 
mic halo, cannot be originated by the complex shift of any 
parent spherical distribution. 
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APPENDIX A: INTEGRAL REPRESENTATION 

We start from the integral representation of the complexified 
potential 



-l>c 



-G 



piv) 



■d?y 



G 



J p{y){h + ih)d^y, 



(Al) 

where r = x — y, and p{y) is the parent density distri- 
bution. As in the real case, general informations about the 
behaviour of the potential can be obtained by expansion in 
special function (such as Legendre polynomials, e.g., Jack- 
son 1999) of the integral kernel. A first possibility is to use 
the expansion in Legendre polynomials for ||r|| > ||a||, as 
given by Lynden-Bell (2004b, eq. [1]) in the case of a point 
charge: 



1 



1 °° 
I'fW /- ^ 



-J'n(M) , 



> Ikl 



where [jl is the cosine of the angle between a and r. Here 
we follow an alternative approach, i.e., first we express in 
integral form the real and imaginary parts of "l>c, and then 
we expand the resulting kernels. We define r = a; — y and 
we consider 



|r — ia|| = llrll — ||a 
so that 



— 2i < a. r >= ve 



V 



= Vi\ 



r\\' — \\a 



2^2 



+ 4 < a,r >2 



and, for v ^ 0, 



sin 6 • 



2 < a.r > 



(A3) 



(A4) 



(A5) 



Note that v vanishes only on the ring determined by the in- 
tersection of the spherical shell 1 1 r 1 1 = 1 1 a 1 1 centered on 1 1 a; 1 1 , 
with the plane P containing the point x and perpendicular 
to a of equation < a,r >— 0; on this ring the phase angle 
e is not defined. We caU D the disc ||rji < ||a||, with r G P. 
We now determine the square root yjve^^^'^^^^^ . All points 
r G P and outside the disc D are mapped to points in the 
(ii,^) plane with sinS = and cosS = 1, so that the norm 
in eq. \h3>\ is made positive by fixing fc = and cutting 
the complex plane along the negative axis, i.e. — vr < < tv. 
Moreover, cos(6/2) > for r ^ D, and the bisection formu- 
lae needed to construct the square root are obtained from 
eq. dm, (EH), JMl) 



(A6) 



9 
cos - = 1 


r- 


+ 


r p — a 2 




21- 


. 9 
sm — — - 






\/2 < a,r > 



vWv+ r , 



(A7) 



The only delicate region is the disc D, where cos 9 = ~1 and 
sin6' — cos{9/2) — 0: the function sin(S/2) is then discon- 
tinuous, changing from -j-l to —1 when crossing the disc in 
the direction of a. After some work we arrive at 



\A^ 



h- 



h = sgn(/i) 



r r — a 



V2v 

\/v~\~\\a\\^~~\ 
V2v 



(A8) 



We now consider the behaviour of eq. HA1|I for ||r|| > \\a\ 
and llrll < ||a||. It follows that for t = ||a||/||r|| < 1 



/« = 



H+{t,fi) 



h = 



sgn{fi) H^{t,fi) 



while for t = ||r||/||a|| < 1 

H-{t,fi) sgn{fi)H+(t,^) 

J^ ^ — n~~n — ' /a = n~~n 



(A9) 



(AlO) 



where 



(A2) 



H+{t,fi) = 



H-{t,,,) = 



,Jl-f^ + ^(1 - i2)2 + 4P^ 

A/2^(l-i2)2+4t2^2 ' 

Y^t2 _ 1 + ^(1 _ ^2)2 + 472^ 
V2^(l-i2)2+4f2^2 



(All) 
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Two interesting formulae can now be obtained for the func- 
tions H+ and H-. In fact, separating eq. (|A2[I in its real 
(even powers) and imaginary (odd powers) parts, and com- 
paring them with eqs. (|A9|I , it follows that 



Jf+(t,/.) = ^(-l)"P2„(M), 



sgn{fi)H-{t,fi) = ^(-l)"P2„+i(^) 



(A12) 



_2n+l 



In other words, H+ and H- are generating functions in 
closed-form for even and odd Legendre polynomials, respec- 
tively. We were unable to find this result in the standard 
literature (e.g., Courant & Hilbert 1953, Morse & Feshbach 
1953, Erdelyi et. al 1955, Abramowitz & Stegun 1972, Grad- 
shteyn & Ryzhik 1980, Arfken & Weber 1995), nor computer 
algebra systems such as Maple® and Mathematica® were 
able to re-sum the series in eq. (IA12|) (while able to confirm 
the expansions). 

With the expressions above eq. (|A1[I becomes 



K($c 



-gE(^ 



Piy) P2nin)^3 



\2n+l 



jr||> 



gE 



(-ir 



p{y)P2n+iiM)\\rf"^'d''y. 



(A13) 

Now, for a density distribution decreasing sufficiently fast, 
the integral over ||a; — j/|| < ||a|| is negligible with respect to 
the other in the far field regions (i.e. | |a;| | ^ oo). In addition, 
as the even Legendre polynomials are bounded by unity, it 
results that the leading term in eq. HA13|I is just the potential 
of the parent system, and this explains why the isodensity 
contours of the shifted models in the external regions are 
almost coincident with those of the parent system. 



APPENDIX B: AN INVERSION PROBLEM 

The procedure is best illustrated as follows. Let assume that 
the (unknown) spherically symmetric real function F{r) (po- 
tential or density, the argument is identical), is the par- 
ent distribution of the investigated axisymmetric system 
f{R,z), so that f{R,z) — 5R(Pc) for some (unknown) value 
of the shift parameter a. The z-shift of amplitude a on the 
radius r gives 

Fc = F(rc) , (Bl) 

and all the considerations of Appendix A hold, where now, 
due to the special orientation of a, r1 — r^—a?'-2iaz = ue'*", 

V = ^(r2_a2)2+ 40^22 , (B2) 



2az 



(B3) 

V V 

By evaluating Tc — y've'^'^, from eqs. (|A6[I - (|A7[I and fol- 
lowing comments, it turns out that on the equatorial plane 

(z^O) 



Fc 



F (V-R2 - a^) , R>a, 

F (iiVffl^ - R^) , R<a, 



(B4) 



where R is the cylindrical radius, and the two signs corre- 
spond to the limits z -^ 0" and z ^ 0^ for _R G D. Equation 
(|B4|) reveals that on the equatorial plane Fc is linked to the 
parent spherical potential by a very simple substitution. In 
particular, the function f{R, 0) for R > a gives the whole 
function F{r) (as the variable \/R^ — a^ spans the entire 
range < r < cxj and -Fc coincides with its real part), while 
f{R,0) for _R < a is the real part of the map of F{r) for 
r £ D. In practice, for an assigned f(R,z), we construct 
the putative parent function F(r) = /(\/r2+a^, 0), with a 
free parameter. We then shift (with a shift parameter b) the 
spherical candidate, and check if its real part can be made 
equal to f{R, z) for some choice of a and b. In case of im- 
possibility, it follows that the original distribution cannot be 
obtained from the complex shift of any spherical system. 

We now apply the considerations above to the MN and 
the Satoh (1980) discs, and to the Binney logarithmic halo. 
Simple algebra shows that for the MN disc, the evaluation of 
eq. ([5} for z — Q and the successive substitution R^ — > r^+a^ 
leads to a spherical candidate F(r) given by the potential 
of a Plummer sphere of scale-lenght a^ -I- (1 -I- s)^. A similar 
result is reached for the Satoh (1980) disc, whose normalized 
potential is 



*(P,z) 



sJr^ + z^ + s{s + 2^/TT^ 



(B5) 



where s is the free parameter controlling the disc flattening 
(Binney & Tremaine 2008). From the substitution R^ = 
r^ -f- a? for jz = we obtain again, as spherical candidate 
F{r), the potential of a Plummer sphere of scale-lenght a^ + 
s{2 + s). Thus, if the MN and Satoh discs are generated by 
the complex shift of a spherical model, this is a Plummer 
sphere. However, the complex shift of a Plummer sphere 
does not generate a MN or a Satoh disc (see CGQ7), and 
this close the investigation. In the case of Binney logarithmic 
halo (|27p . the restriction to the equatorial plane and the 
radial substitution leads to 



Infl 



-r^ + a^) 



F{r) 

but its complexification with shift b reads 



In 



K(Fc; 



\{l + r^+a'^-b'^f+4b^z'' 



(B6) 



(B7) 



an expression coincident with eq. (|27|l only in the trivial 
case 5 = 1 and b — a = 0. We conclude that also the Binney 
logarithmic halo cannot be derived from the complex shift 
of any spherical system. 
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